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Abstract. We examine oscillations as a function of Fermi energy in the capacitance 
of a mesoscopic cavity connected via a single quantum channel to a metallic contact 
and capacitively coupled to a gate. The oscillations depend on the distribution of single 
levels in the cavity, the interaction strength and the transmission probability through 
the quantum channel. We use a Hartree-Fock approach to exclude self-interaction. The 
sample specific capacitance oscillations are in marked contrast to the charge relaxation 
resistance, which together with the capacitance defines the RC-time, and which for spin 
polarized electrons is quantized at half a resistance quantum. Both the capacitance 
oscillations and the quantized charge relaxation resistance are seen in a strikingly clear 
manner in a recent experiment. 
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Figure 1. (color online). Mesoscopic capacitor: A cavity is connected via one lead to 
an electron reservoir at voltage V ac and capacitively coupled to a gate with voltage V g . 
Only one contact formed by a quantum point contact permits carrier exchange. The 
voltage Vqpc controlls the transmission through the quantum point contact. 



1. Capacitance and charge relaxation 

We are interested in the dynamics of a quantum coherent capacitor. We ask: if an 
excess charge has been created due to a sudden change in applied voltage or due 
to a spontaneous fluctuation, how long does it take for this excess charge to relax? 
For a classical capacitor the answer is simply the RC— time. Even for a quantum 
coherent capacitor, the answer can be expressed in terms of a capacitance and a charge 
relaxation resistance. However, the properties of these two transport coefficients can now 
differ dramatically from their classical counterparts. We use the term electrochemical 
capacitance to emphasize that the capacitance now depends on the physical properties 
of the system, and use the term charge relaxation resistance to emphasize that the 
resistance might be very different from what is expected classically. 

The system of interest is shown in figure 1. A small cavity is connected via a 
single lead to a metallic contact and is capacitively coupled to a gate at voltage V g . 
The cavity acts as one plate of a capacitor and the gate as the second plate. The 
opening between the cavity and the metallic contact is formed by a quantum point 
contact (QPC). Of interest is the entire transition from very weak coupling to very 
strong coupling when electrons can transmit the QPC with transmission probability 
1. A structure of this type was analyzed more than 13 years ago, by Biittiker, Pretre 
and Thomas [1] who found a charge relaxation resistance which in the case of a single 
channel of spin polarized electrons is given by half a resistance quantum independent of 
the transmission probability of the channel. Interestingly, a recent experiment by Gabelli 
et al. [2] on a structure similar to that shown in figure 1 measured both the in-phase and 
the out-of-phase admittance and found excellent agreement with the predictions of Ref. 
P]. In particular their experiments confirm the quantization of the charge relaxation 
resistance. 
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Much of mesoscopic physics has centered around a few elementary systems which 
best demonstrate a quantum effect of interest. Examples are isolated rings which were 
shown to support a persistent current [31 El El E] , rings with leads exhibit an Aharonov- 
Bohm effect in conductance [H [8j [10] , QPC's [THUS] and quantum Hall bars [13] teach 
us about conductance quantization and quantum dots [H] reveal in a striking manner 
Coulomb blockade and Kondo physics. Clearly the mesoscopic capacitor is also such 
an elementary system: its major advantage is a complete suppression of dc-transport. 
Unlike in systems with dc-transport where very high frequencies are required to make 
dynamic effects visible on top of the zero-frequency dc-response, here the dynamic effects 
provide the leading order of the system response. 

2. The admittance of a capacitor 

Classically one might want to consider the sample in figure 1 as a series connection of 
a geometrical capacitance C and a resistance R. The capacitance C is the geometric 
capacitance between the cavity and the gate and the resistance R is that of the QPC. 
At low frequencies the conductance G(u) = dI(u)/dV(u) can be expanded in a power 
series to give 



The first term describes an out-of-phase response determined by the geometrical 
capacitance C, whereas the second term is the in-phase, dissipative response determined 
by the resistance R. The quantum coherent capacitor of interest here has an ac- 
conductance which can similarly be expanded in powers of uj but instead of the classical 
quantities C and R is determined by an electrochemical capacitance C M and a charge 
relaxation resistance R q , 



The electrochemical capacitance differs from C whenever the small system is not 
efficient in screening electric fields. The charge relaxation resistance R q reflects the fact 
that carriers injected into the cavity spend a time on this "plate" which is too short to 
equilibrate them [lj [15] and for this reason R q might differ from R. The capacitance C M 
is related to the imaginary part of the AC conductance [H UHl HE] and is often measured 
dynamically [21 [T71 HE] but in theory it can also be obtained by differentiation of a 
thermodynamic (grand-canonical) potential [HI [201 [22] • Here we pursue a dynamic 
approach. 

Biittiker, Thomas and Pretre [1] developed a theory of charge relaxation of quantum 
coherent capacitors based on scattering theory and a self-consistent (Hartree like) 
treatment of interactions. Carriers from the reservoir enter the cavity and are eventually 
reflected. If the contact supports iV transverse channels, the reflection amplitudes can 
be described by an iV x N scattering matrix S. Since all carriers are reflected the eigen 
values of this scattering matrix are of the form expi(p n where <p n is the increment in 
phase which a particle in the n-th eigen channel acquires upon reflection. Both and 




(1) 




(2) 
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R q depend on the derivative of such phases with respect to energy d(p n (E)/dE. The 
derivative is taken at the Fermi energy. When multiplied by h the phase derivatives are 
the Wigner-Smith delay times. Alternatively we can say that 

v n {E) = (l/2 7 r)d(j )n (E)/dE (3) 

is the contribution to the density of states in the cavity of the n-th eigen channel. 
The total density of states in the cavity at energy E is v(E) = (l/2ir)Tr(S^dS/dE) = 
Yjn^n^E). Assuming that the potential inside the cavity can be described by a single 
variable and that Coulomb coupling to the gate is described by a geometrical capacitance 
C, Refs. [U [15] find for the electrochemical capacitance 

_ CeME) 

If the Coulomb interaction is weak C ^> e 2 u(E) the electrochemical capacitance is 
determined by the density of states ~ e 2 v(E), whereas if the energy to charge the 
cavity is large the capacitance is essentially determined by the geometrical capacitance 
~ C. Due to quantum interference the density of states v(E) exhibits mesoscopic 
fluctuations from sample to sample: no two capacitances are exactly equal. For a single 
channel connected to a chaotic cavity these fluctuations are discussed by Gopar, Mello 
and Biittiker [23]. For a cavity coupled to many channel QPC's mesoscopic capacitance 
fluctuations are treated by Brouwer and Biittiker [24J and Brouwer et al. |25j . 
For the charge relaxation resistance Ref. [1] finds, 

h En.U 2 n (E) 

Perhaps the most striking prediction of this approach is that if the connection between 
cavity and reservoir permits transmission of at most one spin polarized quantum channel 
(n =1 in equation ([5])), the charge relaxation resistance at zero temperature is equal to 
half a resistance quantum [TJ [15] 

This should be contrasted with the resistance R oc with T the transmission 
probability, which we would expect if inelastic scattering is sufficiently strong to relax 
the carriers in the cavity toward an equilibrium distribution. The prediction is for a 
resistance that is entirely independent of the scattering properties of the channel! Note 
that the predicted charge relaxation resistance is not h/e 2 , the resistance of a perfect 
two-probe single channel conductor but only half a dc-resistance quantum of a single 
channel. 

The experiment by Gabelli et al. [2] measures both the in and out of phase parts 
of the AC conductance of a mesoscopic RC circuit. The sample similar to figure 1 
with a sub-micron sized cavity is cooled to sub-Kelvin temperatures and subjected to 
an ac voltage in the GHz range. A high magnetic field is applied such that only a 
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single spin polarized edge channel is partially transmitted and reflected at the QPC 
Gabelli et al. [2] develop a model which permits to determine and R q using Eqs. 
( )4ll5l) and compare with experiment. In their model the scattering matrix follows from 
Fabry-Perot like multiple reflections in the cavity. 

The results of this experiment are in good agreement with the theoretical 
predictions of in particular they confirm the universality of the single channel 
charge relaxation resistance over a range of voltages over which the transmission 
probability of the QPC varies substantially. Two effects might give rise to deviations 
from the quantized charge relaxation resistance, equation ([6]): First, as the transmission 
probability T through the contact becomes sufficiently small, the time a carrier spends 
inside the cavity will eventually become long compared to an inelastic time. As a 
consequence we expect realistically that for very small contact transparencies there is a 
cross-over from R q given by equation ([6]) to the "classical" result R oc 3> 4. Second, 
for small transparencies, charge quantization (Coulomb blockade) will become important 
and the random phase approach used by Biittiker et al.pQ might fail. The second 
question has already been addressed in recent work by Nigg, Lopez and Biittiker [26J. 
These authors show that the charge relaxation resistance R q remains quantized even 
in the strong Coulomb blockade limit. To demonstrate this, an approach is developed 
which is based on a Hartree-Fock treatment which explicitly excludes self-interaction. 
For many channel contacts, the resulting expressions for R q , like the expressions derived 
from the Hartree approach, can be formulated in terms of density of states or the time- 
delay of particles. 

It is the purpose of this work to present the results for the electrochemical 
capacitance in the presence of interactions treated on the Hartree-Fock level. Below 
we outline how the density of states are obtained within a Hartree-Fock approach. We 
then present results for the mesoscopic capacitance and the total charge in the dot as 
the transmission through the contact increases from zero to one. 

3. Density of states and charge steps 

For a moment consider the cavity disconnected from the lead. In Hartree-Fock the 
isolated dot is described by the effective single particle Hamiltonian H. In the eigen 
basis {|A}} of the Hamiltonian, we have H = I A) (-^1 where E\ are the Hartree- 

Fock level energies. Coupling the cavity via a single quantum channel to the states in 
the leads with tunnel matrix elements W\ gives rise to a scattering matrix [27] 



S(E) 



1 + iK(E) 
1 - iK{E) 



(7) 



with 



K(E) 



E 



(8) 



x 



E\ — E 



where 



r A = ttW* x w x . 



(9) 
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The function K(E) determines the phase 4>(E) which a carrier acquires upon 
reflection through <ft = 2 axctan(K(E)). The density of states of the cavity v(E) = 
(l/2m)S*(E)dS(E)/dE = (l/ir)daictaji(K(E))/dE is then given by v{E) = Ea^a(^) 
with 

i r A 
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{E x -Ef + (E X -E)E 



(10) 



Note that this expression for the density of states differs from the often used Breit- 
Wigner expression. According to equation (TIDl the density of states at every energy 
depends on all the eigen energies (poles), except exactly at the resonance energy E = E x 
where the density of states is u(E x ) = 1/ttY\. We next specify how we treat interaction. 
In Hartree-Fock the self-energy is [28J 



J HF)mi 



E r 



[ID 



Here E c = e 2 /C is a charging energy. The first term in equation ( ITTj) is the Hartree 
energy determined by the occupied states. The occupation number operator of the 
state I in the dot is rii = d\di where di annihilates a carrier in the cavity. The second 
term in equation (11 ID is of central importance for the description of charge quantization 
effects: the second term excludes self-interaction. Note that the self-energy is a matrix 
with dimension M x M equal to the number of states in the dot. In the following 
we neglect the off-diagonal elements of this matrix. These elements account for weak 
inter-level exchange interactions and for the parameters of interest here affect the results 
only quantitatively (26j. With this approximation the self-consistent Hartree-Fock level 
energies are given by 



E x = e x + E c J2(n»), 



(12) 



where (n M ) = / dEf{E)v ll {E) and e M is the bare (non-interacting) energy of level //. 
Eqs. fl 1 Oil 1 1(1 1 21) are a set of self-consistent equations which we solve iteratively. 



4. Capacitance oscillations 

To proceed, we now assume for simplicity that the non-interacting energy levels are 
equidistant in energy with a spacing Aa = A. We also assume that all levels couple 
with the same strength = 7 to the lead. Furthermore we use a result from Ref. 
[29] which relates the coupling constant and the bare level spacing to the transmission 
probability T through the contact, 

7 = — (2 - T - 2v / T^r) . (13) 

This is a local relation and thus also valid in the Coulomb blockade regime. To complete 
our model we connect the transmission probability through the QPC to the Fermi 
energy, T = l/[l + exp {a(Ep — i?o)}]- This expression is appropriate for a QPC with an 
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E =0.5, M=30, A=1 E =0.5, M=30, A=1 

c c 




Figure 2. (color online) Left: Electrochemical capacitance oscillations (solid blue line) 
as a function of the Fermi energy. The (green) curve following the capacitance peaks 
is the capacitance determined by the escape rate (see text). Right: Charge steps as a 
function of Fermi energy. The transmission probability through the QPC is shown in 
all figures as a (red) dashed curve. 

(inverted) parabolic potential [30]. Here E determines the energy at which transmission 
is 1/2 and a is a constant which depends on the curvatures of the saddle point potential. 

In figure [2] we show the zero temperature electrochemical capacitance C M and the 
total dot charge N for different interaction strengths Ec/A = 0.5 and 1.5 as a function 
of the Fermi energy. In all the graphs the probability for transmission through the QPC 
is shown as a dashed red curve. The total number of levels in the dot is M = 30. 
As expected, for very small transmission, the capacitance has sharp Coulomb peaks at 
energies at which transfer of a carrier into the dot is permitted. The charge increases step 
like with flat plateau regions in which charge is quantized. As the contact is opened the 
peak height decreases and the width of the capacitive Coulomb peaks increases. Even 
when the quantum channel is completely open, the density of states, equation ( fTOl) . 
exhibits weak oscillations. This is similar to a chaotic cavity which even for a completely 
transparent quantum channel exhibits density of states fluctuations and thus capacitance 
fluctuations [23l 1211 125] . The existence of oscillations for T = 1 differs from the model 
used by Gabelli et al. [2j where the cavity was treated as a Fabry-Perot resonator. In 
the model of Ref. [2] the capacitance oscillations are proportional to the amplitude of 
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reflection of the QPC. The model of Gabelli et al. [2] is appropriate for electron motion 
along edge states. Clearly, if the QPC becomes completely transparent to an edge state 
there is no interference, thus no density oscillation and thus no capacitance oscillation. 
In contrast the approach used here is suitable if only a weak magnetic field is applied 
(corresponding to a flux through the cavity of the order of a flux quantum). 

The (green) curves connecting the peaks in figure 2 show the capacitance calculated 
with the density of states v(E F ) = \/ix^(E F ). Deep in the Coulomb blockade regime this 
function coincides with the local maxima of the electrochemical capacitance, however, 
with increasing transparency, as neighboring peaks start to overlap, the maxima of the 
capacitance cease to be determined by a single Hartree-Fock level. 

The graphs on the right hand side of figure 2 show that with increasing transparency 
trough the QPC charge oscillations become ever weaker and for large transmission we 
have a nearly continuous increase of the charge in the cavity. While semi-classical theory 
of the Coulomb blockade leads to charge steps which are infinitely sharp in the zero- 
temperature limit, here the sharpness of the charge steps depends on the transmission 
probability. Other systems which exhibit charge steps which are not infinitely sharp 
even in the zero temperature limit are Cooper pair boxes [31] G2] (where the slope has 
been measured [32| [T7J [18]) but also single electron systems like rings with in line or 
side quantum dots [22J. 

5. Conclusions 

In this work we have examined Coulomb blockade oscillations for mesoscopic quantum 
coherent capacitors. We present an effective single particle approach using scattering 
theory and treating interactions on the level of a Hartree-Fock approach which excludes 
self-interactions. This approach preserves quantum coherence and at the same time 
takes effects of charge quantization into account. For a nearly isolated cavity the 
capacitance exhibits sharp and well seprated Coulomb peaks. As expected, with 
increasing coupling between cavity and reservoir charge quantization effects diminish. 
However, interestingly, we find that even for a completely transparent quantum channel 
capacitance oscillations persist. These oscillations reflect quantum coherent scattering 
in the cavity. 

The results presented here can be experimentally tested. Experiments [2] in high 
magnetic fields have already reported clear signatures of a quantum coherent capacitor. 
Our results should be useful in finding such signatures in the low magnetic field range. 
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